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EXECUTIVE SUMMARY 


This study demonstrates the feasibility of estimating link travel] time and speeds in real- 
time for a rural, mountainous section of an interstate freeway in Colorado using vehicles 
instrumented with global positioning system receivers, that serve as probes in traffic 
streams. The system configuration proposed and developed includes a cost-effective, 
portable GPS deployment unit, communication links to a PC server, and an integrated 
prototype system for vehicle tracking and estimating statistics on link travel time. The 
system also includes a user interface. 


The section of the I-70 corridor in Colorado selected to evaluate the system is located in a 
mountainous terrain just east of the Eisenhower Tunnel, gateway to the ski resorts, 
including Winter Park, Vail, and Aspen, from the Denver International Airport. It carries 
a high traffic volume during the ski season and during the summer months. Currently 
there is very limited surveillance infrastructure to monitor traffic conditions on this 
critical corridor. 


This study, funded by the Colorado Department of Transportation and carried out in two 
phases, explores the feasibility of monstonng this corridor using mobile surveillance 
technology. The objective of the first phase was to identify a cost-effective means of 
monitoring traffic within a rural, mountainous stretch of the I-70 corridor. In addition, an 
algorithm was developed to estimate average link speed, travel time, and the standard 
error of estimates that provides the user with information on the reliability of the 
estimates. The algorithm was tested based on simulated data for different traffic 
conditions and geometric characteristics of links. 


In the second phase of the project, a real-time demonstrational prototype was developed to 
receive, process and estimate link travel time and speed statistics in real-time based on the 
algorithm proposed in the first phase. The system was tested off-line and on-line based on 
field data received from the 1-70 corridor. In addition, GPS receivers were deployed using a 
commercial van pool service. This report presents an overview of the demonstration prototype 
developed, the surveillance technology selected, an evaluation of the system based on 
simulated and field data, procedures to install and maintain the system, and recommendations 
for larger scale implementation. 
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OVERVIEW OF THE REPORT 


The first chapter of this report presents an overview of the specifications of the GPS receivers 
and the modem selected. It also includes a description of the form of communication selected 
to transmit the GPS data from the field to a fixed-end central server. 


The second chapter of the report presents details of the components of the real-time 
demonstration prototype system developed to track probe vehicles in the traffic stream, 
methods to estimate link travel time and speed based on probe data and the graphical interface 
of the system. 


The third chapter presents an overview of the network selected to conduct both on-line and 
off-line tests based on simulated and field data to evaluate the demonstration prototype and the 


findings. 
The fourth chapter presents conclusions and recommendations of the study. 
The appendices include instructions on how to install, run and maintain the prototype system, 


the source code for the programs written, sample GPS data, sample output of the program, and 
the study proposal. 


CHAPTER 1: 
SELECTION OF SURVEILLANCE 
TECHNOLOGY 


L1 GLOBAL POSITIONING SYSTEM (GPS) 


Global Positioning Systems (GPS) allow instantaneous position and velocity of a moving 
vehicle to be determined. To provide continuous global positioning capability, 21 evenly 
spaced satellites placed in a circular, 12-hr orbit inclined at 55° to the equatorial plane provide 
the desired coverage. This constellation provides a minimum of four satellites in good 
geometric position, 24 hours a day, anywhere on the earth. Depending on the selected 
elevation angle, more than the minimum number of satellites are often available for use. 


GPS satellites provide the capability of determining location in terms of longitude, latitude and 
elevation by the simple resection process using the distances measured to satellites. The space 
coordinates relative to the center of the earth of each satellite can be computed from the 
ephemerides broadcast by the satellite using an algorithm. The ground receiver, defined by its 
geocentric position vector, employs a clock that is precisely set to GPS time. The true distance 
range to each satellite may be accurately measured by recording the time required for the 
- satellite signal to reach the receiver. Each range defines a sphere with its center at the satellite. 
Therefore, using this technique, ranges to three satellites are needed since the intersection of 
three spheres yields three unknowns (longitude, latitude, and height). GPS receivers use an 
inexpensive quartz clock set approximately to GPS time. The clock timing error or clock bias 
is overcome by measuring the distances to four satellites. Other errors include satellite position 
error (ephemeris error), ionospheric and tropospheric refraction, receiver noise, multi-path and 
selective availability, collectively referred to as the UERE, or User Equivalent Range Error. 
The cumulative UERE totals are multiplied by a factor of 1 to 6, a factor that represents the 
Dilution of Precision, or DOP. DOP is a measure of the geometry of the visible satellite 


constellation. The ideal orientation of the constellation would require all satellites to be 
equally spaced around the receiver, and one directly above. This would result in a low DOP. 
Selective availability is the highest source of error, an intentional error, imposed to limit 
accuracy to 95% probability of 328 ft. (100 m.) or less. 


Instantaneous velocity of a moving vehicle is determined using the Doppler principle of radio 
signals. Due to the relative motion of the GPS satellites with respect to a moving vehicle, the 
frequency of a signal broadcast by the satellite is shifted when received at the vehicle. This 
measurable shift is proportional to the relative radial velocity. The radial velocity of the 
satellite is known, therefore the radial velocity of the moving vehicle can be deduced from the 
Doppler observable. The accuracy of velocity is + 0.16 kuy/hr (0.1 mph) for receivers with 
location accuracy of 100 m. (328 feet), and with selective availability turned on. Since May, 
2000 selective availability has been turned off and the accuracy has been reported to increase 
tenfold. The instantaneous GPS velocity reporting is independent of position fixes. Prices for 
GPS receivers range from $100-$50,000, depending on the additional features available and 
the reporting accuracy of the receivers. The GPS receiver selected for this study are capable of 
reporting location within a 100m accuracy, providing data in a non-proprietary format and is 
compatible with a Cellular Digital Packet Data (CDPD) modem (Trimble Navigation Limited 
1996). 


Figure 1. Placer 450 and Magnetic Antenna. 


1.2 MODEM 


The Cellular Digital Packet Data (CDPD) is a wireless, public access, packet data standard 
designed to operate over existing analog cellular phone systems. The Uniden Data 1000, a 
CDPD modem, was selected for this study. The selection of a modem is based on: 

« Cost 

» Reliability 

« Compatibility with the GPS receivers selected. 


Figure 2. Uniden Data 1000, CDPD Modem. 


13 COMMUNICATION 


CDPD is a connectionless network service that overlays the cellular voice network to 
wirelessly transmit data to mobile and fixed-end computing devices. The CDPD network 
operates as an extension of the TCP/IP data communication network. TCP/IP is the method by 
which data on the network is divided into packets of bytes. Then the CDPD network transmits 
the GPS data from the cormidor to a PC server in an end-to-end system. 


Key components of the end-to-end system used are the mobile system and the fixed-end host 
system. The mobile equipment, installed in the vehicle, comprises of a GPS receiver and a 
CDPD modem. The fixed-end host system, a PC, runs FleetVision (Trimble Navigation 
Limited 1998b). The PC server is connected to the internet and therefore the PC connects to 


4 


the CDPD network through the internet. Messages received by the cellular ffom the GPS are 
processed and routed to an IP address. Additionally, the cellular carrier providing CDPD 
services accepts TCP/IP transmissions for IP addresses for units that are registered and operate 
in its CDPD network and passes those messages on to the units. The GPS receivers come with 
a SLIP interface driver to communicate with the CDPD modem and a TCP/IP stack to enable 
it to exchange messages over the CDPD network. The CDPD modem is registered with a 
CDPD service provider and configured with the IP address. The modem comes with software 
to allow easy setup. Figure 3 below shows the configuration of the communication network. 


CDPD Modem GPS Receiver 


Figure 3. Communication Network. 


1.4 PORTABILITY 


The GPS receiver, CDPD modem, and their antennas assembled in a box provide a portable 
unit as shown in Figure 4. This umit is easy to install in a vehicle by slipping the adapter into a 
lighter/power source and mounting the magnetic antenna on the rooftop of a vehicle. 


Figure 4, Assembled GPS Unit. 


LS COST 
The cost of the portable GPS receiver unit shown in Figure 4 was less than $1500 per unit 
during the study. Currently, the same unit could be assembled for less. The breakdown of 


costs per unit is presented in Table 1. 


Table 1. Cost Breakdown Per Portable GPS Unit. 


Cost for Study Current Costs 

“GPSReceives (ti (tst(‘<‘w @BWOO”OCOCOCO*O*O*~*~*C*C*«SSSO™ 

CDPD. Modem $350 $350 

Antenna Kit $180 $45 

Cables $60 $50 

Antenna for Modem $50 $50 

Cigarette Lighter Adapter $45 $45 

Cash Box $10 $10 
‘TOTAL — si—‘C;C™~™~*C=«SSSAIS CSA 


CHAPTER 2: 
DESCRIPTION OF THE DEMONSTRATION 
PROTOTYPE 


2.1 COMPONENTS OF THE PROTOTYPE 


The real-time demonstration prototype system developed to recetve GPS data from probe 
vehicles in real-time, process the data and estimate link travel time and speed consists of three 
integrated modules. They are: 


1) Probe Vehicle Tracking Module 
2) Link Travel Time/Speed Algorithm 
3) Graphical Interface 


Figure 5 shows a schematic of the components. A brief description of each module is provided 
below. 


TRAVEL TIME 


¢ Converts GPS time to local time, 
e Estimates ink travel ime and 

| speed over a given interval of 

5m or T5unin 


Figure 5. Components of the Real-Time Demonstrational Prototype System. 


2.2 PROBE VEHICLE TRACKING SYSTEM 


The Probe Vehicle Tracking Module receives GPS data in real-time from all vehicles 
instrumented with GPS receivers or probe vehicles operating within the test corridor using the 
communication network described in the previous section. This information is received as a 
data string at pre-specified time intervals and includes: Report Time and Position Fix Time in 
GPS time, latitude, longitude, altitude, horizontal and vertical speed, number of satellite 
vehicles used, and satellite vehicle ID. The navigational message may be received every 
second. For the prototype, the reporting interval is set to 5 seconds. This data is stored in a 
database file (*.mdb) and is passed to the Travel Time Estimation Module. Probe vehicles 
may be tracked and displayed on a map of the network using Fleet Vision (Trimble Navigation 
Limited 1998a). 


2.2.1 Fleet Vision 


FleetVision (Trimble Navigation Limited 1998b) tracks a fleet of vehicles equipped with GPS 
receivers in real-time. It also reports and stores the data, such as vehicle location and time, to a 
base computer via the CDPD modem. The reported GPS data contain vehicle ID, longitude, 
latitude, altitude, and speed for each vehicle every pre-specified interval. The reporting 
interval can be set to as frequent as 1 second. For this test network, the reporting interval was 
set to 5 seconds. 


2.3 TRAVEL TIME AND SPEED ESTIMATION ALGORITHM 


The link travel time and link speed are calculated by using a program written in Fortran 
language. 


2.3.1 Fortran Program 


The Fortran program has been developed to : 

» Extract the GPS data for a selected network. 

» Convert GPS time to Denver local time. The daylight saving time change and leap 
year are also taken into account. GPS time is the number of seconds elapsed after 
00:00:00 AM January 6, 1980. 

« Estimate the link travel time and link speed based on data received from probe 
vehicles that complete their journey through a link. The link travel time and link 
speed can be determined as follows: 


LTT, =Touy = Tiny (1) 


where, 
LTT, = Travel time for probe vehicle i for link j 


Tig = Time probe vehicle i enters the link j 
Toeg = Time probe vehicle i exits the link j 
L = Length of link j 
tspp, (2) 
LTT, 
where, 


LSPD, = Travel speed for probe vehicle i for link j 


The probe vehicles completing their journey on a link within a given interval are used to 
estimate average link travel time and average link speed. For example, to estimate 15-minute 
average link travel time, all probe vehicles that complete their journey through the link are 
included. Probe vehicles are checked every 2.5 minutes, that is, the update mterval is 2.5 
minutes. 


OT = 3) 


TSPD, =. (4) 
where, 
LTT, = Average travel time for link j 
LSPD, = Average travel speed for link j 
n = number of probe vehicles for link j 


For any given interval, the exact time a probe vehicle enters and exits a link is estimated by 
first converting the longitude, latitude, and altitude at any time to Earth Centered Earth-Fixed 
(ECEF) coordinates and estimating the distance traveled. This procedure is described in the 
next few sections. 


2.3.2 Transformation from Geodetic Coordinates to Earth Centered, 
Earth-fixed (ECEF) Coordinates 


The GPS reports data in longitude and latitude. To calculate the linear distance of the vehicles 
between two report locations, the longitude (¢), latitude (1), and altitude (/) are converted to 
Earth Centered Earth-fixed (ECEF) Coordinates ((x,y,z) coordinates) (Figure 6). 


Figure 6. Geodetic Coordinates Transformation. 
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The flatness of the geodetic ellipsoid is defined as: 


f -— = 0.0034 (5) 
where, 
a = Semi-major axis length = 6378137.0 m 
b = Semi-minor axis length = 6356752.3142 m 
The eccentricity of the ellipsoid is defined as: 


e=/f(2—f) =0.0818 (6) 


The length of the normal to the ellipsoid, from the surface of the ellipsoid to its intersection 


1a Te () 
The linear position in (x,y,z) coordinates is calculated as: 

x, = (N +h,)cos(A,)cos(¢,) (8) 

¥, = (N +h,)cos(A, )sin(¢,) (9) 

2, =(N(Q-e*)+h,Jsin(4,) (10) 


where, 
A, = Latitude at time t 
g,= Longitude at time t 
h,= Altitude at time t 
2.3.3 Estimating the Entry and Exit Time For Probe Vehicles 


Link j-] | Link j | Link j+1 


Figure 7. Probe Vehicle Entering and Exiting a Link. 
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The locations of vehicles in longitude (¢), latitude (2), and altitude (h) are converted to 
(x y,Z) coordinates based on equation (8) , (9) , and (10). The time a probe vehicle enters and 
exits a link is estimated based on linearly interpolating between two consecutive 5-second 
reporting of location as follows: 


h=[(s, 5) +04) Ha -2) Tf cin 
1, =|(%,-%, J +(¥. “Vy, y +(z, -2,) (12) 
1 
b=[le. 2a.) +.) (6. -20.) } (13) 
Ly=|(%,~ 8.) #(0, Pea) +, Aaa) | (14) 
where, 
a = Distance a probe vehicle travels between time t, and t, 
Le = Distance from upstream node for link j to the location of a probe vehicle at 


time tz 
a = Distance a probe vehicle travels between time t, and tpi 
L, = Distance from the location of a probe vehicle at time t, to node }+1 


t, = Last reporting time before entering link j 
t = First reporting time after entering link j 
t, = Last reporting time before exiting link 7 


t = First reporting time after exiting link 7 


A, +424, >%,, = Latitude at time 11,to,ta,and ta, respectively 
¢, %, , $.., = Longitude at time ty, t,t,,and tat, respectively 
Aish, = Latitude for node j and j+1, respectively 


p+ 


9, Pt = Longitude for node j and j+1, respectively 
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Therefore, the time a probe vehicle i enters and exits the link j are as estimated follows: 


Tay 


zs _ yj 15 
4£4+G WT ( ) 


L, 
Teaty = f, eal ( Ar mam (16) 
Based on equation (3), (4), (15) and (16) the average link travel time and speed are estimated. 
In addition, the standard error of the average link travel time estimates are estimated based on 


the average variance and covariance of the probe vehicle estimates as follows: 
2 J 
oT (4,)=¥; (4)+ yl (t,)-¥, (1,)] (17) 


Or (¢,) = standard error of average travel time estimate for link j, for a given analysis 
period (e.g. S-min or 15-min) 

v,(t,) = average variance of average travel time for link j, , for a given analysis 
period 

n,(t,) = average covariance of average travel time for link j, , for a given analysis 
period 

n(t,) | = number of probe vehicles that travel link j, for a given analysis period 


Similarly, the standard error for average travel speed estimate is also calculated, All estimates 
for every link m the network are stored in a file (*.csv) and passed to Graphical Interface 
Module. 


2.4 GRAPHICAL INTERFACE 
To display the average link travel time and speed estimates, their statistics and individual 
probe vehicle data in a user-friendly environment, a graphical interface was developed. The 


routines of this module use the Open Database Connectivity (ODBC) interface to receive data 
from the Travel Time Estimation Module. 
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2.4.1 FleetPlot 


FleetPlot has been developed as an Excel application. The FleetPlot is used to plot vehicle link 
travel time, vehicle link speed, vehicle instantaneous speed, average link travel time and 
average link speed. The plots show the data for the test network. The plots are automatically 
updated every 2.5 minutes. This update can be stopped or restarted. The plots only show the 
time of the last update and the number of probe vehicles in a given analysis period. The plots 
or data display are delayed approximately by a maximum of 2.5 minutes. Figure 8 shows a 
snapshot of the Graphical Interface. 
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Figure 8. Sample Interface Showing the Integration of the Three Modules. 


25 PROTOTYPE INTERFACE 


The prototype consists of three modules: Probe Vehicle Tracking, Travel Time Algorithm, 
and Graphical Interface. The three modules are connected using Open Database Connectivity 
(ODBC) interface. The Travel Time Algorithm module receives data from the Probe Vehicle 


Tracking system. The Graphical Interface module receives data from the Travel Time 


14 


Algorithm as shown in Figure 5. The external program performs all the calculations such as 


time conversion, link travel time, link speed, etc. (Figure 9). 


Automatically update every 2.5 minutes 


Probe Vehicle Tracking Database 
(avil_00.mdb) 


Graphical Interface module imports GPS data from Probe Vehicle 
Tracking database by using ODBC interface (Open Database 

Connectivity) and exports data as CSV files. The CSV Giles serve as 
mput for Travel Time Algorithm module. 
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Travel Time Algonthm module performs the following tasks: 
© Time conversion 

© Calculate link speed et 

@ Calculate link travel time [ 

@ Calaniate average link speed and link travel time for a 

vehicles within {5 mins. interval, | 

® Process output files Le. linst.csv, [link.csv,.. 4 


ane Oe ce ny nny Gamer Om ne ey 


Graphical Interface module imports the data by using ODBC ss 
interface and plots graphs. 
® Travel time for link j for probe vehicle no. 1 to 2 
© Speed for link j for probe vehicle no. 1 ton 
@ Vehicle speed for vehicle no.] ton 
© Average travel time 
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Figure 9, Flowchart. 
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CHAPTER 3: 
EVALUATION OF THE DEMONSTRATION 


PROTOTYPE 


To perform a comprehensive evaluation of the demonstration prototype, a section of the I-70 
corridor between the Denver International Airport and the ski resorts was selected. The next 
few sections present details of the test section selected, a simulated representation of the test 
network, off-line tests conducted, on-line real-time field data collection procedures, and 


findings of these tests. 


3.1 TEST NETWORK 


Interstate 70 westbound between US 40 Exit (Winter Park) and US 6 Exit was selected as the 
test section. The length of the 2-lane test network is approximately 15.8 miles. The average 
grade is about 3% upgrade, with 2 maximum of 6%. The on-ramps and off-ramps within the 
test network include US 40, Georgetown, Silver Plume, Bakerville, Herman Gulch, and US 6. 
Figure 10 shows the test network. Figure 11 shows the elevation of the network. 
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Figure 10. Test Network. 
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Figure 11. Elevation Profile. 
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3.22 NETWORK REPRESENTED IN SIMULATION 


To develop and test the algorithms under different flow conditions, probe vehicle deployment, 
and probe report interval, the test network was simulated using CORSIM, a microsimulation 
model. Field data collected may often be limited. Therefore, a simulation provides extensive 
data for testing under a variety of operating conditions. The simulation model was calibrated 
based on field data to ensure that the network was well represented by the simulation model 


3.2.1 CORSIM 


CORSIM (CORidor SlMulation) (Federal Highway Administration 1996), a micro- 
simulation, developed by the Federal Highway Administration (FHWA) consists of NETSIM 
and FRESIM and is used to simulate a surface street and freeway network, respectively. 


The traffic environment of the network is represented to include the topology of the roadway 
system, roadway geometry, lane channelization, driver behavior, .traffic control, traffic 
volumes entering the roadway system, turn movements, and fleet characteristics. A schematic 
of the representation is shown m Figure 12. 


3.2.2 Calibrating CORSIM 


To calibrate CORSIM, traffic flow was recorded at US 40 and Bakerville using video cameras 
as shown in Figure 13. A video image processing system was used to estimate traffic flow rate 
and speed at these locations. In addition, volunteer drivers with GPS receivers were deployed 
to collect travel time data along the same section as flow data was recorded. 
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Figure 12. Nodes and Links Represented in CORSIM. 
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Figure 13, Cameras at US 40 and Bakerville to Monitor Traffic Flow on the I-70 Corridor to Calibrate CORSIM. 


19 


3.2.3 Results of Calibration 
Details of the calibration results are presented in several papers. ((Khan (in preparation)), 
(Khan 1998), (Khan 2000b), and (Khan 2000a)). 


3.3. OFF-LINE TESTING OF THE AVERAGE LINK TRAVEL TIME 
ESTIMATION ALGORITHM 


Analysis of the standard error of the estimates of average link travel time based on Eq. (17) 
and probe vehicle data shows that increasing the number of probes beyond a certain level does 
not improve the accuracy of the estimates as shown in Figure 14. Here the curve flattens as the 
number of probes traveling a link increases, with very little marginal improvement in 
accuracy. In addition, the marginal improvement in the estimate of average link travel time 
also varies based on flow conditions and link characteristics. Details of these findings are 
reported separately (Khan 2000b). Figure 14 shows the standard error for a 833.63 m. 
(2,735 ft) link with a 3% grade between Bakerville and Herman Gulch for two flow 
conditions. Figure 15 shows the standard error for a short weaving section between an on and 
off-ramp. Figure 14 and Figure 15 show that the marginal improvement in error is higher for 
the weaving section for less than 5 probes, compared to a longer section. 
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Figure 14. Standard Error of Average Link Travel Time Estimate for Link between 
Bakerville and Herman Gulch. 
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Figure 15. Standard Error of Average Link Travel Time Estimate for 
Link between US 6 On and Off-Ramp. 


Therefore, using link travel time functions, average link travel time estimates based on probe 
data for an analysts period, and the number of probes used for the estimate, the standard error 
of the average travel time estimate may be reported. For example, for probe data collected 
within a 15-minute analysis period and using Eqs. (1) to (16), the average travel time for Link 
# 316-416 may be estimated. Based on the link travel time function as shown in Figure 16, the 
flow level may be identified as “Medium”. Using standard error curves such as Figure 15 and 
the number of probes within the analysis period, the standard error of the link travel time may 
be estimated. 
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Figure 16. A Link Travel Time Function. 


3.4 OFF-LINE AND ON-LINE TESTING OF THE DEMONSTRATION 
PROTOTYPE SYSTEM BASED ON CONTROLLED EXPERIMENTS 


Several vehicles were equipped with the assembled GPS unit (Figure 4) and volunteer drivers 
were deployed for several weekends: March 25, Apnil 1, and April 15, 2000. The data 
collection was performed during 7:00AM to 10:00AM. During this period traffic volume 
varied from 700 to 2200 vph westbound for the test network. 


For the field study, the test network was divided into 10 links based upon the location of on- 
ramps and off-ramps. The length of links varies from 0.12 to 4.12 miles. The upstream and 
downstream locations and lengths are as follows: 


US 40 off-ramp to US 40 on-ramp, 1.11 mites. 

US 40 on-ramp to Georgetown off-ramp, 3.16 miles. 
Georgetown off-ramp to Georgetown on-ramp, 0.22 miles. 
Georgetown on-ramp to Silver Plume off-ramp, 2.11 miles. 
Silver Plume off-ramp to Silver Plime on-ramp, 0.16 miles. 
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Silver Plume on-ramp to Bakerville off-ramp, 4.12 miles. 
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7. Bakerville off-ramp to Bakerville on-ramp, 0.21 miles. 

8. Bakerville on-ramp to Herman Gulch off-ramp, 2.73 miles. 

9. Herman Gulch off-ramp to Herman Gulch on-ramp, 0.24 miles. 
10. Herman Gulch on-ramp to US 6 off-ramp, 1.73 miles. 


3.4.1 GPS and Modem Coverage 


Extensive testing of the prototype system bas shown that the GPS coverage varies on the test 
network on I-70 depending on the season, with better coverage during the summer months. 
Based on a 5-second reporting interval, the GPS coverage statistics are summarized in 
Table 2. As shown in this table, the average consecutive seconds either the GPS signal or 
modem coverage is lost due to inadequate satellite coverage, mountainous terram, or any 
obstructions is less than one minute during winter and approximately half a minute during 
summer. The average total duration signal is lost for a one-way, 16-mile trip is a little over two 
minutes during winter and less than two minutes during summer. 


Table 2. Statistics on Loss of Probe Vehicle Data for the Test Network. 


Average Average (o 
consecutive iON Si r duration signal 
seconds GPS or or modem 


modem coverage | connection for a | modem coverage | connection for a 


is unavailable is unavailable | one-way trip for 
(seconds) test network is (seconds) fest network is 
lost (seconds) lost (seconds) 
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3.5 ON-LINE, REAL-TIME TESTING BASED ON COMMERCIAL 
VANPOOL SERVICE AS PROBES 


There are several van pool services from Denver International Airport (DIA) to ski resorts 
such as Colorado Mountain Express (CME) (Figure 17). For this study, Colorado Mountain 
Express was contacted to explore the feasibility of deploying several assembled GPS units 
from DIA to Vail to test the demonstration prototype in real-time. CME has 150 vans in their 
fleet operating on this route with an average headway varying from 20-30 minutes to an hour 
for winter and summer season respectively. CME agreed to participate in this study and 5 
units were deployed between June 27, 2000 to July 17, 2000. 


Figure 17. CME Vanpool Service. 


Extensive on-line testing revealed a few software problems that were later resolved as part of 
this process. The demonstration prototype is now debugged and is fully operational and 
running in the Colorado TransLab. 


CHAPTER 4: 
CONCLUSIONS AND RECOMMENDATIONS 


The study demonstrates the feasibility of estimating link travel time and speeds in real-time for 
a rural, mountainous section of an interstate freeway in Colorado using vehicles instrumented 
with global positioning system receivers, serving as probes in traffic streams. The system 
configuration proposed includes a cost-effective, portable GPS deployment unit, 
communication Jinks to a PC server and an integrated prototype system for vehicle tracking, 
estimating statistics on link travel time and displaying the information in a user friendly 


environment. 


Extensive testing of the prototype system has shown that the GPS coverage varies on the test 
network on I-70 depending on the season, with better coverage during the summer months. 
Based on a S-second reporting interval, the GPS coverage statistics show that the average 
consecutive seconds either the GPS signal or modem coverage is lost due to inadequate 
satellite coverage, mountainous terrain, or any obstructions is less than one minute during 
winter and approximately half a minute during summer. The average total duration signal is 
lost for a one-way, 16-mile trip is a little over two minutes during winter and less than two 


minutes during summer. 


The study also shows that average link speed, a measure independent of link length, is a better 
estimate to display and for an operator of the system to quickly draw inferences about traffic 
congestion. Link travel time estimates are also required for targeted traveler information for a 
given origin and destination pair. The design of the most efficient and user-friendly front-end 
to the prototype was not part of the scope of the study. This study may be further extended to 
address these issues. 
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A final project presentation was made on July 20, 2000 to the Project Manager, traffic 
engineers from the Mobility Group, DTD, ITS Group, Headquarters, Region 1 and Region 6, 
CDOT, and the Denver Regional Council of Governments (DRCOG). The prototype system 
presented was very well received by the attendees. Interest was expressed by several groups 
within CDOT to further extend this study for other applications. The groups that expressed an 
interest included: GIS, Mobility, ITS, Region 6 and Region 1. Several follow-up presentations 
were scheduled. One of the meetings arranged by the Mobility Analysis Group on August 
23 2000 included 20 attendees from the GIS, Mobility Analysis, ITS, and Region 6, CDOT. 
A separate meeting was scheduled by Region 1 to be held in October, 2000. 


42 EXTENSION OF THIS STUDY 


A proposal was submitted in December, 1999 to CDOT as part of the 2001 Problem 
Statements to extend this study to use the Regional Transit District (RTD) bus fleet already 
equipped with Automatic Vehicle Location System (AVL) to estimate freeway traffic speed. 
The proposal was very well received and has been recommended for funding. This new study 
is titled "Using RTD's Transit Vehicles as Probes to Develop Speed Maps for Colorado 
Freeways". Bruce Coltharp from the ITS Office, CDOT will serve as project manager. 
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APPENDIX A: 
INSTALLING AND RUNNING THE PROTOTYPE SYSTEM 


B.1 


SYSTEM REQUIREMENTS 


« Pentium 200 MHz class or higher with CD-ROM drive 
» 96 MB of memory 

e 2.0 GB hard disk drive 

s Microsoft Windows 98 

» Microsoft Excel 97 

« A fixed-end host server with a static IP address 


B.2 INSTALLING FLEETVISION 


B.2.1 


» Insert the FleetVision installation CD in the computer's CD-ROM drive. 
= Click on Start and select Setting/Control Panel. 

» Click Add/Remove Programs. 

» Inthe Add/Remove Programs Properties dialog box click on Install and 


follow the on-screen instructions. 


Communication Channel Setup 


1. In FleetVision Taskbar, Click Tools and then Channel Setup. 
2. Data layer should be set to TAIP as shown in Figure A-1.. 


Figure A-1. Channel Setup Dialog-Data. 


3. Protocol layer should be set to TCP/IP CDPD and the IP address is for the fixed- 
end host that will receive the GPS data. 


Figure A-2. Channel Setup Dialog-Protocal. 


B.2.2 Vehicle Configuration 
1. In FleetVision taskbar, Click Fleet, Vehicle Configuration, and then Create 


a new Vehicle or edit an existing vehicle. 
2. Click Continue... 
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Figure A-3. Configuration Manager Dialog. 


3. Enter Physical ID (a number to represent the vehicle), Name and select icon by 
clicking Icon. 
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4. Enter communication information by clicking Setup to access to Vehicle 
Communication Configuration Dialog. 
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Figure A-5, Vehicle Communication Configuration Dialog. 


5. Select CDPD in Select Protocol and then click Setup to enter the CDPD 
modem JP address. 
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B.2.3 Installing FleetPlot and Fortran Program 


1. Insert the FleetPlot&Fortran program installation CD in the computer's CD- 
ROM drive. 

2. Click on Start and select Setting/Control Panel. 

3. Click Add/Remove Programs 

4. Inthe Add/Remove Programs Properties dialog box click on Install and 
follow the on-screen instructions. 


B2.3.1 MS-Dos Prompt 


MS-Dos prompt must be set to minimized as follows: 


1. Click Start, Programs, and then MS-DOS prompt. 
2. Click £7] MS-DOS Prompt Properties dialog. 
3. Set Run to Mmimied. 


| MS5-D05 Prompt 
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Figure A-6. MS-DOS Dialog. 
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Figure A-7. MS-DOS Prompt Properties. 


B.2.4 Installing GPS Receivers 


B241 Installing GPS in a Probe Vehicle: 


1. Place GPS assembled unit (Figure 4) the GPS box in your car close to the 
cigarette lighter socket. The GPS unit may be placed on the front seat. The 
data modem's antenna should be in an upright position. . 


On-Off Power Switch 


Figure A-8. Data Modem. 


2. Attach the magnetic GPS antenna on probe vehicle's roof 
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Figure A-9. GPS Unit (Placer 450) and GPS Antenna. 


3. Phuig the power cable to probe vehicle's cigarette lighter socket. 
4, Turn the modem power switch on. 
5. Check the data modem by following these steps: 
« Power LED- Solid red indicates power on and no light indicates 
power off. 
+ TX LED indicates modem operational status. 
e Solid Orange indicates trying to acquire a CDPD channel. 
» Blinking Orange indicates trying to register on an acquired 
CDPD channel. 
¢ Solid Green indicates registered. 
«Flashing Green indicates transmitting while registered. 


While operating, check the CDPD modem to make sure: 
Power LED - is red all the time. - 
TX LED - ‘agree mos of thetime, 


B33 SETTING UP THE MODEM 
The fixed-end host computer recetving the probe vehicle data must be on the internet service 


with a static IP address. The power of the modem should be on. In order to set up the modem 
to report GPS data to the computer, the following steps should be performed: 
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Click Start, and then click Run. 

Type Telnet and then click OK. 

Click Connect and then click Remote system. 

Enter the modem IP address. e.g. 166.130.7.3 

Enter port 59100 and Terminal type vt100. 

Click Connect 

Type > sdaOcOxxx. xxX.xxx.xxx< , replacing xxx.xxx.xxx.xxx with the fixed - 

end host computer's IP address 1.e. > sda0c0132.194.014.007 < 

8. Type >qda< and then press Enter. A message should be displayed similar to 
the following. >RDA0C0132.194.014.007;ID=0003< 

9. Close telnet dialog box. 
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Figure A-10. Teinet Dialog Box. 
B.4. RUNNING THE PROTOTYPE 


To nun the prototype, the Fleet Vision and FleetPlot program should be open while the external 
program is called by FleetPlot. The windows of two programs: FleetVision and FleetPlot will 
be displayed. The FleetVision tracks the vehicles. The FleetPlot displays the plots. 


B.4.1 FleetPlot 


In order to start FleetPlot, the macros should be allowed to run by clicking Enable Macros. 


Figure A-i1. Excel's Macros Protection Dialog. 


B411 = FleetPlot menu 


The charts show the vehicles traveling on I-70 from exit 232 (Winter Park exit) to US 6 off- 
ramp for a test network. The number of probe vehicles is currently set to 5 vehicles. 
Additional probe vehicles may be added as needed. 
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Figure A-12. Fleetplot Menu. 


Link Travel Time and Link Speed 


Under the FleetPlot menu, this selection allows a user to view the bar charts of link travel time 
(minutes) and link speed (mile per hour) for probe vehicles 1 to 5. The link travel time and 
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link speed charts show data for the last westbound run. For any link if the data from the last 
run is not available, the data from previous run will be shown instead. The charts also show 
the last updated time for each link or the last reported time before the vehicles move across the 
link to the next link. The link travel time and link speed for each link is represented in different 


color. 
Vehicle speed 


Under the FleetPlot menn, this selection allows a user to show a plot of speed of the probe 
vehicles every 5 seconds. Each color will represent the data for each link. There are currently 
10 links on the test network. The charts show the GPS speed data for the last run on WB I-70. 
The upstream links are on the nght of the chart and downstream links are on the left of the 
chart. 


Average Last 15 min. 
Under the FleetPlot menu, this selection allows a user to show the average speed and average 


travel time of the vehicles for each link withm a 15-minute analysis period. The number of 
vehicles on the link within a 15-minute period is also shown on the table below the chart. 


StopAutoUpdate 


Once the FleePlot starts, the data is automatically updated every 2.5 minutes, This function 
may be disabled by clicking "StopAutoUpdate" 


RestartAutoUpdate 


Under the FleetPlot menu, this selection allows a user to enable the auto update function The 
data may be automatically updated every 2.5 minutes. 


Update Now 

Under the FleetPlot menu, this selection allows a user to start update function immediately. 
The updated plots are displayed after 2.5 minutes. The automatic update function is disabled. 
To automatically update the charts, click "RestartAutoUpdate”. 


Network Map 


Under the FleetPlot menu, this selection allows a user to show a map of the test network on 
I-70 from US 40 off-ramp to US 6 off-ramp. 


Help 


Under the FleetPlot menu, this selection allows a user to show information about the FleetP lot. 


A-10 


APPENDIX B: 
MAINTAINING THE PROTOYPE SYSTEM 


C.l1 FLEETVISION AND FLEETPLOT MAINTENANCE 


The Probe Vehicle Tracking module runs FleetVision (Trimble Navigation Limited 1998a) 
and its database. The Graphical Interface module runs a program FleetPlot. Depending on the 
frequency of probe data collection, the database may require regular maintenance, The 
maintenance routine requires the probe data to be exported as text files and saved. The 
schedule for maintenance depends on the number of probe vehicles in the fleet. The following 
steps are recommended: 


« Click Fleet on Fleet Vision taskbar (Figure B-1). 
« Select Database Management and Export or Delete data from the database and 
then click Continue button in the Configuration Manager dialog box.(Figure B-2). 


2S ; os 
Figure B-2. Configuration Manager Dialog Box. 
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C.1.1 Export Data 


«s To export data, click 4 | button (Export Data to Text Files) in the Database 
Administrator dialog box (Figure B-3). 

» Enter the export date for one day of data in the Export Date (Figure B-4) and 
then click OK. 


Figure B-4. Export Date Dialog. 


C.1.2 Delete Data 


« To delete data, on the Database Administrator dialog box, enter number of days 
i.e. enter O days to delete all the previous data from database. 


§ button on Database Administrator dialog box to delete data. 
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« When finished, Click Done on Database administrator dialog box. 


C.1.3 Compact AVL Database 


2 Click Start, and then click Programs. 
» Click FleetVision, and then click Compact AVL database. 


To maintain the database, perform the following on a regular schedule: 
¢ Export data from database to text files. 

« Delete data from database. 

« Compact AVL database 
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APPENDIX C: 
EXTERNAL PROGRAM 


| Lastchange T 13Dec1999 1:44pm 

C Lastchange T 6Sep2000 11:25am 
PROGRAM LINK TRAVEL TIME 

C This progmm is used to perfixm the fhllowing make: 

c = Extract WB data on the test network 

c¢ — -Canvert GPS time to local time. Daylight saving hes also 

¢ been taken into account. 

¢  ~Caloulate Link travel time, link speed, their avemee, and instaniancous 

¢ speed 
CHARACTER IMOIDJREPTIME,LATILONGIJALTISPEED 
1 JNFILE(5)*12, OUTFILE1()"9, OUTFILE2(5)"9,OUTFILEX5)°9 
1, TEMP1(5)"1L, TEMP2(S)*1 1, TEMP(5)*10 
CHARACTER*25 REAL TIME(9999), LIRTIME($999), LIRTIME(9999) 
1,L3RTIME(9999),L AR TIME(9999),L SR TIME(0999) LER TIME(999) 
1,L7RTIME(9999), LER TIME(9999), LOR TIMIE(9999), LIOR TIME(9999) 
1,221,222,273, 724 295 796 277 ZR TIO) 
REAL*8 MOI(20000) REP(20000),LA(20000),0(20000),ALT(20000) 
1,SPE(20000),rep1(20000) 
1, TMO1(18000), TREP(18000),TLA( 18000), TLO(1 8000), TAL T(18000) 
1, TSPE(18000),MOID1(9999), REPTIME 1(9999),LAT1(9995) 
1, LONGI (9999), ALTITUDE 10999),SPEED 1(9999), RLONG( I) RLAT(11) 
1,TD(1),PDC).TTNG TAD 
LLTTILSPI,LTT2,LSP2,LTTALSP3LTTALSP4LTTSLSPSLTTELSPS 
LLTT7LSP7.LTT&LSPRLTT9LSP9LTTIO,LSPIO 
1,A,BX1,Y1,C1 P1P2,P3 PP LPP? PP3 lo(20000) dt 1(20000) 
1,DLO2(20000) 
1,LISPD(10800),LZSPD(1 0800),L3SFD(10800) 
1, L4SPD(10800),LSSPD(1 0800), LASPD( 10800) 
1,L7SPD(10800),L8SPD(I 0800), L9SPD(10800),L 10SPD(10800) 
1 DIFTIME diffe 
112073 2AZ5, 6278 I9,210 
1,211,Z12,213,.Z14Z1SZ16 217 Z18219,.20 
1,231,Z30,733,734, 235, 26,737,238 739740 
1, HDEN(1.0800) HIDEN72(10800) HDEN'(1 0800) 
1,GPSHOUR(10800),GPSTIME( 10800), DENH( 10800) 
1 XCD) RVC DARZ(1)X(1G800), ¥(10800),2(10800) 
1,DLAT1(.0800), DLONG) (10800), DHEZGHT(10800), TN(20800) 
INTEGER DENYEAR(10800), DENMONTH(10800), DENDATE(1 0800) 
1, DENHOUR(10800), DENMINUTE(10800) DENSECOND\( 10800) 
1. MOTH(13),J, NN(10) 


c 
c 
C PARTI EXTRACTING WEST BOUND LAST RUN DATA 
c 
Cc 


DOLOOP 5 TIMES DUE TO THERE ARE 5 INPUT FILES(S VEHICLES) 


DO NNNN=1,5 
C FILEi={NPUTFILE 
C FILE 10 AND 1 ARE OUTPUT FILES, THE REST ARE TEMP. FILES 
INFILE()}= "VEHICLEI.CSV” 
INFILE(2)="VEHICLE2.CSV" 
INFILE@)= "VEHICLE3CSV" 
INFILE(4= “VEHICLEACSV" 
INFILE(5= “VEHICLESCSV" 
OUTFILEI(1)="UINST.CSV" 
OUTFILE1(Q2) = "2NST.CSV" 
OUTFILE}(@3)="SINST.CSV" 


C-} 


OUTFILE}(4)="4INST.CSV" 
OUTFILE1(5)="SINST.CSV" 
OUTFILEX1)=*ILINK.CSV 
OUTFILE2(2)="2LINK.CSV” 
OUTFILE23)="SLINK.CSV" 
OUTFILE2(4)="4LINK.CSV" 
OUTPILEXS)="SLINK.CSV" 


C OUTFILESS ARE THE FILES THAT CONTAIN LINK SPEED, LINK TRAVEL 


TIME,AND 


C CLOCK TIME FROM PREVIOUS UPDATE. THE DATA FROM PREVIOUS 


UPDATE WILL BE USED IF 
C THE LINKS CONTAIN PARTIAL DATA (POR PRESENT RUN, LINK 
TRAVEL TIME AND 


C LINK SPEED WILL BE CALCULATED FOR ONLY FULL DATA LINKS. 


TEMPI(1)="IWBLAST.CSV” 
TEMPI(2)="2WBLAST.CSV" 
TEMPI(3)="SWBLAST.CSV" 
TEMP1(4)="4WBLAST.CSV" 
TEMP1(5)="5WBLAST.CSV” 
TEMPA1)="IWBALLR.CSV" 


TEMP2(5)="SWHALLR.CSV” 
TEMPX(!)="IWBLRT.TXT” 
TEMPX2)="2WELRTTXT" 
TEMP33)="3WBLRT.TXT” 
TEMPX4)="4WBLRT-TXT" 
TEMPX5)="SWBLRT. TXT” 
OPEN(FILE~INFILE(NNNN)) 
OPEN( FILE = TEMP! (NNNN)) 
OPENGFILE=TEMPUNNNN)) 
OPEN(GS FILE=TEMPANNNN)) 
OPEN(OFILE= OUTPILE(NNNN)) 
OPEN(1,FILE=OUTFILEX(NNNN)) 
OPEN (12,FILE = OUTFILEXNNNN)) 
C BELOW ARE LONGLAT OP EACH NODE 1 =US400FFRAMP,10= 
US6OFFRAMP 
RLONG(I)=-105647973 
RLONG(2) =-105663047 
RLONG(3)=-105696491 
RLONG(4)= -105696940 
RLONG(5)=-105724720 
RLONG(6)=-105727350 
RLONG(7) = -105803200 
RLONG(B)=-105806990 
RLONG(9) =-105855212 
RLONG(10)=-1058 59005 
RLONG(11)=-105885278 
RLAT(I) =39761585 
RLAT() =39751714 


RLATG) = 39717156 
RLAT@) =29714134 
RLATY5) =39605444 
RLATI(6) = 39654400 
RLAT() = 39693105 
RLATIS) = 39692285 
RLAT@) =39701951 
RLAT(10)= 39701809 
RLATIN 1) =39685803 
C (%\YZEARTH CENTERED EAR TH-FIXED(ECEF) COORDINATES FOR 
NODE1TO 11 
RX{1) -1324786.3 7064 104:R (14729559253 16RI6 
R2(1)-4058779, LOWTI IG 
RYQ2- 1326224 2795ZER YO)--4 729900428 
R27}=4088449. 37683902 
RAGT-1329662.1 1125545 RVI F-4731535 96093 
RZ{3)-4055540,1 1206678 
RX(4)-1329756 46763008 R Y(4)-473 1729.3 S04904 
R2(4)4055278.15985 105 
BX(5)~1332482.455307BGR YU5)~-47325 1 t 6525 160 
R25 405381 t 28702252 
RX{6)=-1332688,85773879;R Y(6)=-47325 18.50449557 
R2(6)=40537 19.2851 9985 
RX(T)+ 1539019. S466 ER YT ATIOGE) 2698 79) 
RZ(T}=4053 736 SBOERIA 
RX(8y~1339349.52256452:R Y(B)~4 730957. 79373887 
R2Z(B)\=4OS36H.4SHOTSAS 
RX(9)~ 10175 BADDSIT7ZE V(9)—47 292203 $86 7H 
R2O)=4054594.1 7085953 
RX(10)>-1343550.26335670.RY (1074729193, 70795796 
R2Y10)=4054587,0 1094854 
RX 1)-1349993, 4559220, RY (1 1 4729797 290 
R211)-4053285 58720716 
C ENDOF REFERENCES LONGLAT 
C RESET ALL THE VALUES TO ZERO 
DO I= 19999 
MOKD-OREPOOLAM)-OLOT), ALTI-O.SPT)-O 
TMORD-O-TREPYIF0, TLAQ)=0, TLO( 0: TAL TO, TSPEXIR=O 
MOIQ)-O.REP(PGLAM CLO ALTO GSPEG)R0 
LISPOG)=0.G.L25PD()-0.0,L ISPD) =0.0,LASED (10, O;LSSPDO-0,0 
‘L6SPD1)+0.0,L78FD0)=0 OL 8SPD9)=0.0,L SSPD1}0.0,L10SPDG0.0 
LIRTIME()>™L2R DMEG)>"L3R TME()-" LAR TIMED 
LSRTIME()="L@R TMEQ="*{L 7 TIME()-™ Lak TM 
LORTIME()="*;L1R TIME()"* MOD 100.0, reptinel ()-0 
END DO 


REAIX1*)IMODDREVTIMEILAT ILONGJALT SPEED 
Cc PABTU . 
C INPUTDATA WERE FILTERED TO OBTAIN 
C LONGLATIN STUDY AREA, AND TO DISTINGUISH WB FROM 
C WBEB 
MO 
kab 
DOI= t, 20000 
READ(L*, END=24MOIM REP().LAM.LOMALTO. PED 
Line below will be weed if PosFix Timo ja comstent(Vehiche 3 
oot moving), the number of seconds aqual to report aterval will 
be sdded to PosFbcime. report aiterva] change, this nted to 
be changed as well. Currently, set to 5 soponds. 
if Gp) .2q. rep(l-1)) then 
rep Q)-epid-1+3.0 
fee 


ae fF 4 © 
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repl@=rp® 
aif 
rep) = rept} 
dint) =0 
dolg=0 
IF (LOO) LE. RLONG(1) AND. LO() G2. RLONG(! 1) THEN 
IF (LAM) LB RLAT() AND. LAQ) GE. BLAT(1)) THEN 
DloG) = (LO0)-L09-1) 
401) = LO0)-L003)) 
iad{]) = (LOG-3}LOW-#) 
= diffime=sep ep 1-1) 
F (dbo) LT.0.0) AND, dio}() LT.0.0 AND. DLO2().LT.0)THEN 
WRITES 25)MOM)REPIDLAM. LOM ALT) SPE 
M=Mil1 
20END IF 
ENDIF 
END F 
24 ENDDO 
25 FORMAT (6(F143,",) 
C FILB3 CONTAINS WADATA FOR ALL RUN 
C BELOW IS TO OBTAIN 70-WB ONLYLAST TRP 
NEWRUN=1 
REWIND (3) 
DOT=14M 
READG,*, END=27) IMOH(), TREPQ), LA), LOG), TALI) TSPED, 
TREP(O)=0 
DIFTIME = TREPQ)-TREP(-)) 
& diftime below is the time betwee rams IN ORDER TO GET ONLY LAST RUN 


I Q@IFTIME GT. 1500 AND. TLO) LE TLO(I-1)) THEN 
NEWRUN =I 
XNDP 
JE (TLOQ) G8. TLOG-1)) THEN 
NEWRUN=1 
END F 
MOID1(D= THOM) 
REPTIME)()- TREPO) 
LATIQ}= TLAQ) 
LONGI()= TLO@) 
ALTITUDEI(= TALT) 
SFEED (Ip TSPE() 
Z7ENDDO 
REWIND) 
c 
NI-O.N2-0,N3=0:,N4-0.N5=O,.NE-O: N70, NS=NO=ON 10-0 
NN(1) NN(Z)=O:NNGYONN(4)20,NN(S)PONNGGH0 
NN(7) =G,NN(8)-QNNO)0,NN(10)°0 
DOJ=NEWRUNM  - 
c¢ GROUPING THE DATA TO EACH LINK(WITHIN THE BLOCK) 
C UNK 1 US40GFFE-RAMP AND US40 ON-RAMP 
IF (LONGI()-RLONG(2).GE. 0) _AND(LONG 1Q)-RLONG(1) .LE.O )jthen 
WRITE(2,25MOID10) REPTIME]()LATIO).LONGIG) 
tL ALTITUDE! (SPEED 1(J) 
NISNI+1 
END F 
C UINK 2US40 ON-RAMP AND GEORGETOWN OFF-RAMP 
IF (LONG1Q) LE. RLONGQ)) AND{LATIG).GE. RLAT(3)) THEN 
WRITE(2,5)MOID1(@), REPTIME)(), LATI(), LONGIG) 
1 ALTITUDE) (),SPEEDI() 
N2=NOQ+] 
END IF 
C LINK 3 GEORGETOWN OFFP-RAMP AND ON-RAMP 


IF ((LATIQ) LE RLAT[3)) AND. (LATI() GE RLAT(4))) THEN 
WRITH(2.252401D} (1), REPTIMEI(), LATIA), LONGIQ) 
1 ALTITUDE) SPEED 1() 
NENG 
ENDIF 
C_ LINK 4GEORGETOWN ON-RAMP AND SILVER PLUME OFF-RAMP 
IF (LAT1@) LE. RLATI4)) AND.(LONG1G).GE. RLONG(5))) THEN 
WRITE(2,25$vi0 1D (7), REPTIMIEIZ), LATIA), LONGI), 
1 ALTITUDEI@),SPEED1(1) 
N4doN+l 
END PF 
C CINK 5 SILVER PLUME OFF-RAMP AND ON-RAMP 
IF QLONGII) LERLONG(S)) AND. (LONGI() GE RLONG(6))) THEN 
WRITE(2,25)M@01D1 (0), REPTIMEIQ), LATIN, LONGI) 
1 ALTITUDE] (SPEED 1{1) 
NSNS+1 
END IF 
C UINK 6SILVER PLUME ON-RAMP AND BAKER VILLE OFF-RAMP 
TF (LONGI) LE RLONG()) AND.(LONG}() .G&. RLONG(7))) THEN 
WEITE(2,25)MOID} (1), REPTIME!(), LATIN, LONG(1) 
1 AL TITUS (SPEED! (1) 
NGNGH 
END F 
C LINK 7BAKERVILLE OFF-RAMP AND ON-RAMP 
IF (LONGI) LE RLONG{7)) AND.(LONG!() GE RLONG())) THEN 
WRITE(2,25)MOID1 (0), REPTIME?(J), LAT), LONGI) 
| ALTITUDE! ()SPRED1() 
NEN 
END F 
C LINK 8 BAKERVILLE GN-RAMP AND HERMAN GLUCH OFF-RAMP 
IF ((LONGIQ) LE RLONG({8)) .AND.(LONG1(@) GE. RLONG@))) THEN 
WRITH(Z,25)MOID 10), REPTMEI@), LATI(), LONGI0) 
1 ALTITUDS1(),SPEEDI() 
NS2NSt] 
ENDIF 
c LINK 9 HERMAN GLUCH OFF-RAMP AND ON-RAMP 
IF QLONGIQ) LE RLONG(9)) AND(LONGI().GE RLONG(I{)) THEN 
WRITE(2.29)MOID! (1), REPTIMEIG), LATICO), LONGI) 
1 ALTITUDBI), SPEED 1(1) 
NON] 
END F 
C LINK 10 HERMAN GLUCH ON-RAMP AND US6 
IF (LONGIQ) LE RLONG(10)) AND(LONGI() GE. RLONG(11))) THEN 
WRITEQ,259)MOID 10), REPTIMEXD), LAT1(T), LONGI) 
LALTITUDE]@) SPEED IC) 
NION1O+1 
ENDIF 
K =MAX(N1N2N3B.N4NSNGNTNBNONIO) 
NNQ)=NI 
NNQ@) = NNUL)EN2 
NNQ) = NNQ2}B 
NN(4)~NNG)HNA 
NING) <NN(4)ENS 
NN(6) = NN(SPANG 
NN(7}=NN(GHNT 
NNG) = NNC7ENB 
NN(S) = NN(BY+ND 
NN(10)=NNOPNIO 
C NNI=NUMBER OF POINTS AT END OF LINK 1. NN2-END OF LINK? 
ENDDO 


REWIND) 
C ABOVE=> WB DATA FOR LASTRUN WBRE SAVED IN FILE 2 
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C eee 
C PART It CONVERT GPSTIME TO DENVER TIME 
Cc 


IF (NN((0) BQ. 0) THEN 
GOTO 200 

END F 

DOI= 1, NN(10) 

READ(2,$, END=30)MOKD,REP(D.LA)LOMAL TD, SPAM) 

ENDDO 


C &OF SECONDS TO 01/01/97 1200;00AM+7ERS.=536137200 SBC OR 


148927HRS, 
¢ Day ight saving: between Ist sunday of April and last sunday of Ocenber 
¢ add mothe | hour 
30 DOI=), NN(10) 
GPSTIME()=REP() 
GPSHOUR()) =GPSTIME(Y3600- 148927 


GOTO 2 
END F 
BDENT)=DENHQ 
END DO 
32 DENYEAR()= 19MIN 
HDEN2() = GPSHOUR()-HDEN() 
C D@NYBAR DEFINES YEAR 
C HORND =4 OF HOURS(MAX =# OF OLR IN THAT YRAR) 
DOIK=3.MN4 
IF (K 2Q MN) THEN 
M=) 
ELSE 
M~0 
END 
END DO 
C IFM®@1JNDICATES LEAP YEARG66 DAYS A YEAR) 
if (MEQ 0) THEN 
MOTH()=00 
MOTRO)~ 744.0 
MOTAQ) = 1416.0 
MOTHG)= 2160.0 
MOTH(4)= 2880.0 
MOTHS) = 3624.0 
MOTH(6)~ 4344.0 
MOTEKT)= 5088.0 
MOTHS) 5832.0 
MOTH(S) = 6552.0 
MOTHQ0)~ 7296.0 
MOTH(11)=8016.0 
MOTH(12)= 8760.0 
ENDIF 
IF (4 EQ. 1) THEN 


MOTH(O)= Q0 
MOTH()— 7440 
MOTH) = 1440.0 
MOTAQ) = 2184.0 
MOTHEA) = 2904.0 
MOTHS) = 3648.0 
MOTH(S) = 4368.0 
MOTEY?)= 5112.0 
MOTH(S) = 5856.0 
MOTHO) + 6576.0 
MOTH(0)= 7320.0 
MOTH()1)= 8040.0 
MOTH 2) = 8784.0 
END IF 
DOL=112 
DENMONTH(D=Z 
IF (ADEN2() LT. MOTH(L)) THEN 
GOTO 14 
END IF 
END DO 
C HDEN3>4 OF HOURS IN THAT MONTH(MAX=#0F DAY IN THAT 
MONTH? 24) 
34 N=DENMONTHO}1 
ADENG() = HDENZ)-MOTHIN) 
DENDA TRQ) = INTEHDENS(1y724)+1 
DENHOUR() = INTGIDENX1)+24-DENDA TE()°24) 
DENMINUTEQ) =INT((HDEN3()-DEND A TR(1)? 24424) DENHOUR(D)*60) 
DENSBCONDO)=NINT((SIDENS()-DENDA TE(}*24+24}- DENBOUR())*60 
*DENMINUTE())"60) 
Add | hour duc to daylight saving 
if Genmonth(i} ge. 4 and. denmomibf) Je. 10} then 
denbour(i) = denboun()+1 
endif 
if (denyean()) 27, 2000) then 
if (denmowth{i) eq. 4 nd. dendateGi) tt 2) then 
denhowst) — desabou(} i 
endif 
if denmonth§) 0g, 10 2nd. dendatofi) gt. 29)then 
denthour())* denhourf)-1 
endif 
endif 


A 


if (Geryoo() .eq. 2001) then 
if (dervennth(}) .og. 10 and. dendate(i) pf 28) then 
denhour() = denhonr(i)-1 
end if 
oud if 


if (Gemyean()) 29, 2002) then 
if @enmouth() .29. 4 and. denate()) fb 7)then 
deaboui) = denhan()-1 
endif 
iF (deneocusiirff) 9. 10 red. dendaw)) gt 27) then 
derheun()= denbaur()-] 
end if 
end if 


if (demyearG) 2. 2003) then 
if (demmanth{)) oq. 4 2nd dendatofi) t.6) than 
denhour@) = denhaur()-1 
end if 
if (densnoethf) oq. 10 and. dendato(f) pe 26) then 


CA 


if (denyenn() eq. 2004) then 
if (denmouth() eq. 4 20d. dendate() 4) then 
denhourf) = denhar(i}-1 
end if 
end if 


if (denyusr(i) .e9. 2005) then 
if (denmanthG) .og. 4 and. denclesefl) it, 3) then 
denbomr) ~ denbourG)-) 
end iT 
if Genatenth(i) 29. 10 and. dendate}) gt 30) then 
_ denbour()+ denhour(i-1 
ond if 
end if 


if (denyern() -<x). 2006) then 

if Genmondy()).eq. 4 and. dendutg) br 2) then 
denhcan()= denbow6)-) 

ond if 

if Genmionth() eq. 10 and. dendatefi) gt 29) then 
denhowG)= denhauG)-) 

endif 

end if 


if (denyear{) 2g. 2007) then 
if (denmnonth() .eq. 10 and. dendate(}) gt. 28) then 
denhourG) = denbour(i)-] 
endif 
endif 
c 
if (GearyearG) .eq, 2008) then 
if (denmonth(?) .0g, 4 snd. dencete() ft 6) then 
denhomG) = denbourt)-1 
endif 
if (denmanth()) 29. 10 and. doodate() gt. 26) then 
denbour(?) = denhour6)-? 
entit 
endif 
c 
if (denycas(t) 2g. 2005) then 
if (Germmcamb()) .2q.4 and. dendate{) in S) then 


if (denmoath() og. 10.and. dencaie(?) gt. 25) then 
denboun(i) = denbourG)) 

end 

endif 
Cc 
if (donyoest) £9, 2010) then 

if (denenenth() .2q. 4 -and. dendate(?) 5} thea 
denhour() = denhouriy! 

endif 


35 FORMA NTLF5.,TIOP}4.3, T26F143,T41F143,TSQFIOZTIOFIOI 
1, THLL 7. T8312,"/, T8646, T9 LI2,"*,T942*?, T9712) 
ENDF 
IF (DENMONTHQ) GE: 10) THEN 
WRITE(5,36)MOKD,REP(),LA(), LOM), ALT(D,SPEQ) DENMONTH() 
1. DENDA TE), DEN YEAR(), DENHOLR() 
LDENMINUTE(), DENSECOND() 

6 FORMARTLES OTIOF 4.3, T25,F 143, TSF 14.3, TSF 10.2, TALE 103 
1,T8L22,"7, T8427, T8714, 19222"2,195.2,"2,T98D) 
END F 
ENDDO 


C ENDOF TIME CONVERTING:WB LASTRUN DATA WITH DENVER TIME 


WERE SEVED IN FILE $ 


¢ 
C PARTIX COMPUTELINK TRAVEL TIME AND LINK SPEED 


c 
C NODE 1=U840 AND NODE 11= US6 
REWIND (5) 
DO 1= 1.NN(10) 
READ (5,49 END=S0)MOID 1 ()) REPIIMEL(),LATHILONG1() 
LALTITUDE]@), SPEED 1()REAL TIME() 
49 FORMAT(TLFS.0,TIOF 14.3, T25F 143,741 F 143, TS6F 10.2, TOF 103 
1,78),A25) 
ENDDO 
C CONVERTINGLONG, LAT TO X Y,Z(METERS) 
DOI= LNNOO) 
DLATI(= LATIQ)"2. 14 1592654/ 1801000000 
DLONGI()= LONG 1(f)*3.141592654/180/1000000 
DHEIGHT(D = ALTITUDEI(1V100.0 
DX) =6378137.04SOR T(I-0.08 18° 2*SINQDLATI(D)**2))) 
XQ) © (IN) +DHEIGHT() *COSMLATICY*COSDLONGIM) 
YQ) = (CIN@D +DAEBIGETOYCOSDLATI(D)SIN(OLONG] (1) 
Zt) ~ CNG) *(-0.0818°*2)-DEEIGHT() *SINDLA TIO) 
END DO 
¢ CONVERTING INSTANTANEOUS SPEED FROM CM/S TO MPH 
JP (NN()) NE. 0) THEN 
50 DD 1= LNN() 
LIS?0q) = SPEED 1)*3600'160934.4 
LIRTIME() = REAL TIMEQ) 
END DO 
ENDIF 
IF (NNQ@) NE. 0) THEN 
DOT=NN(1)+?, NNC2) 
L2SPD() = SPEED 1{1)"36001 60934.4 
LIRTIMAD) = REALTIME() 
END DO 
ENDIF 
1F (NNG) NE. 0) THEN 
DO l= NNG@}LNNG) 
L3SPD(D = SPEED L(D36001 60984,4 
LARTME() = REALTIME() 
ENDDO 
ENDIF 
IF (NN(4) NE. 0) THEN 
DOI= NNG)+1.NNG4) 
LASPD()) = SPEED K(1)*3S00/1 609344 
LARTIME() = REALTIME() 
END DO 
END 
IF (NN(5) NE. 0) THEN 
DO I= NNG@HLNNG) 
LSSPD(])= SPHED 1(1)"360/ 1609344 
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LSRTIME()=REALTIME() 
END DO 
END IF 
IF (NN@6) NE 0) THEN 
DO I> NNO} INN(G 
L6SPDG) » SPEED 1(1)*3600/1 605044 
LER TIME) = REALTIME() 
END DO 
END IF 
IF (NN().NE. 0} THEN 
DOT=NNG@}H1, NN(7) 
LISPD(1)= SPEED 1(1)°3600/160034.4 
LR TIME() =REALTIME() 
END DO 
END F 
If (NNQ) NE. 0) THEN 
DO1=NN(7}+ 1, NN) 
L8SPD() = SPEED 1(1)*3600/180934.4 
L&R DMEM) = REALTIME(D) 
END DO 
END IF 
IF (NN(9) NE. 0) THEN 
DOI= NNG@}HI, NNO 
LOSPDA1) = SPEED 1(0)*3600/160934,4 
LOR TIMB()=REALTIME() 
END DO 
END IF 
IF (RN(10) NE. 0) THEN 
DO I=NN(@)}+1, NN(I0) 
LIOSPD()= SPEEDI()*3600/160934.4 
LIOR HIME() = REALTIMEQ} 
END DO 
END F 


C OUIPUTINITANTANEOUS SPEED FOR EACH LINK 


WRITE(10,51)"MOID","NO. 1", “LINK I TIME", "US40(OFF)-US40(0N)" 
1,°NO.2";"LINK2TIMB’,"US40(ON)-GEORGE TOWNIOFFY" 
1,°NO.3°"LINK3TIME”, “GEORGE TOWN(OFF)GEROGE TOWNUN)” 
1,"NO.4","LINK4TIME*, "GRORGE TOWN{ON)- SILVER PLLUME/OFFY 
1JNO.2"LINKSTIME™ “SILVER PLUME(OFF)SILVER PLLUME(ON)" 
1,."NO.6,"LINK6TIME",“SILVER PLUME(ON)- BAKERVILLE(OEFY" 
1,°NO,7°,"LINK7ITME*,"BAKER VILLE(OFF)-BAKERVILLE(ONY” 
1,"NO.8°,"LINKSTIME", "BAKER VILLE(ON)-HERMAN GLUCH(OFF)" 
1,"NO.9“LINKSTIME”, “HERMAN GLUCH(OFF}HERMAN GLUGHKON)" 
L™NO.10", “LINK IOTIME", “HERMAN GLUCH(ON IL B(OFF)" 
L*US4QQOFF- RAMPS 

$1 FORMAT 31(A°,").A) 
DOIF 1K 
T=NN{1y4 
R=NNQY}I 
D=NNG}I 
YW=NNGH 
G=NNGH 
b~NNOH 
T7=NN()4A 
B=NNGHI 
D=NNOYH 
110 =NN(IO}H 
WRITE. 10,52)MOID (1).LL IRTIME(D,LISPDQ) 
1,31,L2R TIMEG1),L2SPD¢1) 
1,2.LIRTIME(Z)L3SPDL2) 
1,,LAR TIMBQ3)LASPDG3) 
1,4, LSR TIME(4)LSSPD(H) 


1S5,LER TIME(S) LESPOXIS) 
1S. LIRTIMEMO)I SPD6) 
17, LAR TIME) L8SPO(17) 
1 J8,LOR TIMES), L9SPLX13) 
LYLIOR TIME \L 10SPD(9) 
S2FORMATEF143,10((010),,A “AMPM?,.P 14.3) 
ENDDO 
C CALCULATE TIME FOR BACH NODE BY INTERPOLATION 
C TD2=TOTAL DISTANCE BECGPS LOCATION BEFORE AND AFTER 
NODE=ONODE2 
C PD2<DISTANCE BET. GPS BEFORE NODE? AND NODE2 
C TIN2=TIME BRIWEEN GPS BEFORE NODE2 AND NODE2 
Pl=0 
Po 
P=0 
C NODE] AND NODE 1) ARE CALCULATED BY BXTRAPOLA TION 
TOG SQRTCK XUAN VEY 242-21") 
A =(X(NN(10))-X(NN(I)-1)}°°2 
B= (Y{NN(10)}-XQNNGO}-1))**2 
C= ZRRN(IOQ-ZONN(I O71) )°*2 
‘TOI SQRT(ATBHC) 
IF LONG KNN((D) EQ. 0) THEN 
TX =0 
END F 
C PDt=DISTANCE FROM NODE 1(US400FFRAMP) TO FIRSTPOINT 
C TD) =DmSTANCE BEIWENN 1ST AND 2ND POINTS OF GPS DATA 
PDO FSQRTWRAI)-XMYe24RY (IYO)? *2HRAT-Z1)°*2) 
PP) RKC 1)X(NN(1O) )P°2 
Pr2=(RY(1))-YINN(IO)))°"°2 
PRS = (RZ(11}-2(NNCIO))) 2 
PDQ -SQRIPPHPP2+PP3) 
TF(LONGI(NN(0)) EQ, 0) THEN 
POI) =0 
ENDIF 
TTNG)AREPTIME}(2) REPTIME](1)°PDC YTD!) 
TIN(5 1) = (REPTIMENNN(10))- REPTIME](NN{10)-1))"PD(1 1) TEX] 1) 
IF (10(11) EQ. 0) THEN 
TIN() <0 
ENDIF 
C ENDEXTRAPOLATION OF NODES 1 AND Ii 
DOI=19 
IF (NN(H1) NE NN@) THEN 
Cc ** IF THATLINK CONTAINS FULL DATA => LINK TT AND LINK SPEED 
WILL RE CALCULATED 
X11 =(X(NN(D+1)-XQNNG@))*"2 
Y1>(YONNO 1 EYONNID)PS2 
1 = NNO} 1)-ZANNGD)P"2 
TDG 1)-SQRTEI+Y +01) 
PL =(AXGH+1-X(NND))2 
P2=RYCH1)-YNND)))"°2 
Pd = (RZ(H1)-22NN@)Y°2 
PDU+1-SQRTPHP2P3) 
IF LONGINN()) EQ 0.OR LATI(NN() BQ. 0) THEN 
PDH) >0 
ENDIF 
TING 1) = (REPTIME NNN(}+1)- REP TIME L(NNG))) "PD 1 TOG 1) 
IF @DG+)) EQ. 0.0R TH+) EQ. 9) THEN 
TINE 1-0 
ENDIF 
END F 


ENDDO 


READ (1257)Z1,22,21,74 25, 26,Z7,28.29,210 
READ (12,57)Z11,212,713,214,Z15,Z16Z17,Z18,Z19.220 
READ (12,58)221,222,273, 234 225,226,272) 228,229.20 
READ (12,59)231,232,233,234735,236,237,238,239,ZA0 
57 FORMAT (IQF73)) 
58 FORMAT (10(A 19) 
59 FORMAT (IO 14.3)) 
REWIND(12) 
CLOSE(12) 
C —==END INTERPOLATION LINK CONTAINS PARTIAL DATA -LINKTT 
AND LINK SPEED 
C—-——WILL BE SET TO ZERO 
DOI= 
TQ) = REP TIME I(NNG-1)}+ TENG) 
END DO 
IP (PDC) GE. 1000) THEN 
WH)=0 
ELSE 
TU) = REPDMEI(C}-TIN(1) 
ENDIF 


C IF THEREISNODATA ON THATLINK TiN)= NN-1), LITAND LSP=0 
C LTT=LINK TRAVBLE TIME(MINUTES), LSP=LINK SPEED(MPH) 


TP (NU) .NE.0) THEN 

IF (NN!) NEO AND. NNQ).GT. NN(I)JTHEN 
if (tin(1) .ne. 0 pd. tt{Z) ve. 0) then 
IF (1{1) NE. 0} THEN 

LIT (1(2}T11) 60 

LSP 1={S843 A28/S2B0YLTT1°60 
ENDIF 

END if 

ELSE 

Lm =2Z1 

LSP]=Z1l 

LIRTIME(NN()) = 221 
REPTIMENNN(1))=Z31 

ENDIF 

END 


IF (NN(2) NE. 0) THEN 

IF (NNYZ) GT. NNUDDAND. NG) GT. NN(2) THEN 
if (m2) ne. 0 and. tin(3) 10 0) then 

LTT2= (1G)-T2)y0 
LSP2{16725. 5O/SQROVE TI 2" 60 
ENDif 
BLS 

LTT#2 

UsP2=Z12 

LIRTRME(QNN())= 272 
REPTIME1(NNQ@))= 232 
END 

END 


TF (NNX3) NE 0) THEN 
IF (NNG).GTNNQ)-AND. NN(4) GT. NN(G)) THEN 
if (in(3) ve. O and. tin(4) zne. 9) then 
LTT=(14)-TO)VOD 
LSP3~(1139.125/5280)LT13760 
END if 
ELSE 
LTRB=233 
LSPI=Z13 
L3RTIMENNG)) =Z23 


REPTIME)(NNG))~ 233 
END FF 
END F 


IF (NN(4) NE. 0} TEEN 
IF (NN(4) GT. NNG) AND. NN(S).GT. NN{4)) THEN 
if (en(4) te, © and. nS) 24 0) then 
LT TS (1{5)- Tia) y60 
LSP4(L 1158. 39G/S280LTT 4°40 
END it 


END F 
END F 


IF (NNG) NE. 0) THEN 
FF (NNG) GT, NN(4).AND. NN(6) GT. NN(G)) THEN 
dtn(5S) ne. 0 znd. tin(6) ne. 0) then 
LTTS= (TOTS) VO 
LSPS=(838.457/ S2BOYLTTS*GO 
ENDif 
ELSE 
LTIS=Z5 
LSP5=Z 15 
LSRTIMELNNGS)) = ZES 
REPTIMELNN(5))= 235 
END IF 
END IF 


IF (NN() NE. 0) THEN 
IF (NN(6) GT. NN(5) AND. NN(7) GT. NN(6)) THEN 
if (i5x(6) ne. 0 and ten(7) 200. 0) then 
LTTé=(17}- NOY 
LSPO(Z1T76ZIHSISOYLTTSMGO 
ENDif 
ELSE 
LT16=26 
LSPS =Z16 
LARTIME(NN(O)) ~ 226 
REPTIMEI(NNK6)) = 236 
END IF 
END F 


IF (NN() NE 0) THEN 
FF NWNGJ GT. NN(6) AND. NN(8) GT. NN(7)) THEN 
if (tn(7) ne 0 2nd. tl) 18. 0) then 
LTT7= (TR)-IHYEO 
LSPP(1134,640/S280VLTT?°60 
END 
ELSE 
LTT2Z? 
LSPPEZI7 
LIRTIME(NN(T)) 227 
REPTIMEI(NN()}= 237 
END IF 
ENDIF 


FF (NN@®) NE. 0) THEN 
FF NN(S) GT. NN(DLAND, NNO).GT. NN(§)) THEN 
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if (t10(8) ne. 0 and. wm{9) ne, 0) then 
LITS= (TOP Tyan 
LSP8={1442 ),616/S280VLTTE*60 
ENDF 


LSPS=Z18 
LERTIME(NN®)) ~Z28 
REPTIME](NN@)) = 238 
END F 
END F 


IF (NNO) NE. 0) THEN 
IF (NNYS) GT. NN{(8) AND. NNO 6) GT. NN(9)) THEN 
if (an(9) ne. 0 and. an(10) .0¢. 0) than 
LIT9= (TOO}-TY9) 190 
L&P9=(1257.546/5280LTT9°60 
ENDif 
ELSE 
LID-2a 
LSP9Z19 
LORTIME(NN(S) = 229 
REPTIMEI(NN(8))= 239 
END IF 
EOF 


IF QNN(10) NB, 0) THEN 
IF (NN(I0) GT. NN{S) AND. NN(9) GT: NN(B)) THEN 
if (en(16) 00. 0 anid. ttn{11} 00.0) then 
LITNIOTNIIETOeO 
LSP10=(9156.00 ZB0YLTTIO°GO 
END f 
ELSE 
LTTtO>Z10 
LSPIO-Z20 
LIGRTIMEKNN(10))=230 
REPTIMEINN(10)) = Z40 
END IF 
END IF 


6 YF LINK TRAVEL TIME FROM PREVIOUS RUN EQUALS TO ZERO THEN 


C SETITEQUAL TODATA FROM PREVIOUS RUN. 
IF (LTT) BQ, 0.AND. NN(1).NE. 0) THEN 
LT =2) 

LSPI=Zi1 

LIRTIME(NN(1)) = 221 

ENDIF 

IF (LTT2.BQ, 6.AND.NN@Q) NE. 0) TAEN 
LTEZ=Z2 

LSP25Z12 

LIRTMBINN(2)) = Z22 

END F 

YFLTT3 BQ. 0.AND. NN(1) NE. 0) THEN 
LITERS 

LSP3=213 

IR TIME(NN(3))=Z2 

ENDIF 

IF (LYT4 EQ, 0 AND. NN(I) NE. 0) THEN 
LIT4™7ZA 

LSP 214 

LARTIME{NN(4)) = Z24 

ENDIF 

IF (LTTS EQ. 0.AND. NN(1) NE. 0) THEN 


LTTS =ZS 
LSPS-Z 15 
LSRTIMEINNS)} = 225 
END IF 
IF (LTTS EQ. 04ND. NN(1) NE 0) THEN 
LYI6 =2Z6 
LSP6=Z15 
LERTIME(NN(6})) ~Z26 
END F 
IF (LTT7.Q. 0 AND, NNO) NE. 0) THEN 
LTT7427 
LSP7=217 
LIRTIME(NN(7) = 227 
END F 
IF (LTTS EQ.0.AND. NN(1) NE 0) THEN 
LTT =28 
LSPE=Z18 
L3RTIMEINN(S)) =Z28 
END IF 
IF (LTT9 EQ. 0 AND, NN(1) NE. 0) THEN 
LTT9=z 
LSP9=Z19 
LOR TIME(NN(9))= 229 
END F 
IF (LTT10 EQ. 0_-AND. NN(1) NE. 0) THEN 
LITI0=Z10 
LSPle-2a0 
LIOR TMBWNN(16}) =230 
ENDIF 
C SesReeQU PUT Es oesere 
C NODE 1= US400FF NODE |) = US6OFF 
60 FORMATO 1€15.3,..") 
70 FORMAT (10(°,"F15.3)) 
WRITE(1 1,80)" TIME US40-USS,"USAQ(OFEHUS4OQQQNY 
t,"US4Q(0N)-GRG TWN(OFF)” 
1, ,GRG TWN(OFF}-GRG TWNION)* 
1*GRIGE TWN(ON}SLVR PLM(OFP)" 
1,"SLVR PLM(OFF)-SLVR PLM{ON)" 
L°SLVR PLM(ON)-BKRVLL(OFFY 
L"BKRVLL{OFF)-BKRVLL{ON)” 
1,"BXRVLL(CN}URMN GLCHOFFY 
1,"HRMN GLCH(OFP)-HRMN GLCH(ON)" 
1,"HRMN GLCH(ON)-USS{(OFFY" 
SOFORMA RA, 10",",A)) 
90 FORMAT (9(A25," AMPM",","),A25,"AMPM™) 
DOI= 110 
IF (NN) EQ. 0) THEN 
NNO) =1 
LIRTIME() = 0/0/0000 00:00:00" 
L2RTIME() = 00000 00:00:00 
L3R TIMED) = ‘00/0000 00:00:00 
LARTIME() = 010000 00:00:00 
LSR TIMED) = ‘OXy0000 00:00:07 
LOR TIME() = 00000 00,0000 
L7R TIMED) = ‘01/0000 00:00;00 
LARTIME() = (0/0000 00:00:00 
LOR TIME() = OR0000 00200700 
LIOR TIME() = '00/0000 00:00:07 
REPDMEI(}-0 
END F 
END DO 
WRITE, OLTTLLITZLTBLITALTISLTIGLIT7LTVS LTIRLITIO 
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WRITE(11,70)LSP1,LSP2,LSP3,LSP4,LSP5,LSP4,LSP7,LSP8LSP9,LSP 10 
WRITRY LYA,"AMPRPYL 1B TIMEQNN( | J} 
WRITE(13,(A,"“AMPMY)LZR TME(NN(2)) 
WRITE(1L,(A,"AMPM*))L3R DME{NNG)) 
WRITE(11,(A."AMPM"))LAR TIME(NN(4)) 
WRITEILYA“AMPM”)JLSR DME(NN(5)) 
WRITE(I1,(A,"AMPMY'LGR DIME(NN(S)) 
WRITBUEL(A,"AMPM"))}L7R DME(NN)) 
WRITE(11,YA, “AMPM"YILER TIME(NNG)) 
WRITEL1L{A,°AMPM"Y)LIR TIME(NNG)) 
WRITE(ILYA,"AMPMY\L OR TIMEQNN( 10) 
open (12, file= outfile3{nnnn)) 
WRITE(ASTLINLTIALTTALIVALTISLITOLTTLIVALITALIMNIO 
WRITE(12,57)LSPLLSP2.LSPR LSPA LSP 5 LSPRLSP7LSP8LSPSLSP 10 
WRITE( 12,58). 18 TIMBINNG())) LZR TIMENNG)) 
1 L3RTIMENNG)) LARTIME(NN (4) ESR TIME(NN(S)) 
1 LERTIME(INN(6)) LIRTIME(NN(7)) LER TIME{NNGS)) 
LLIR TIMEQNNGD)) L 1GR TMEINN(10)) 
WRITE(12,59)REPTIME I(NN(1)) REPTIME}(NN(2)) 
1 REPTIMBINNG)) REPTIME I(NN(4)) REP TIMEL(NN(S)) 
L.REPTIME}(NN(6)) REPTIME (NN(7)),REPIIME(NN()) 
1 REPTIME 1(NN(9)) REPTIMEI(NN(10)) 
98 FORMAT (XFT.3,"""}F73) 
96 FORMAT (XAI9,,"LAI9) 
97 FORMAT(SXF143.".),F143) 
100 FORMA T(XA 18),A 18) 
CLOSE(Z,STATUS delete’ ) 
CLOSE(3, STATUS = ‘delete! ) 
CLOSE(S STATUS = Welete' ) 
REWIND(1) 
REWINIX}0) 
REWINIXI) 
REWIND(32) 
CLOSE(1) 
CLOSE(10) 
CLOSE(11) 
closel12) 
& chwe(15) 
200 end do 
CALL SMSO 
C ABOVEEND DO IS FOR THE LOOP TO READ 5 INPUTFILES VEHICLE!- 
5.CSV 
END PROGRAM 
¢ SUBROUTINE SMS BELOW IS TO CALCULATE SPACE MEAN SPEED IN 
CASE MORE THAN } VERICLE EXIST 
C ON THATLINK WITHIN 15 MINUTES. 
SUBROUTINE SMSQ) 
DOUBLE PRECISION TIME(S, 10) SPEXS,10) GDIME(S, 10) DD(10) 
1 ATTIME(10), ASPEED(I0). LAS TRUN(10), TTEME(10),N( 10) tspci(20) 
CHARACTER RTIME{5,10)* 19, TR TIMEY) 0" 19 
OPEN (20, FILE > "PILNK ext") 
READ (20, | LOYOME(, 1) TIMB(),2), TIME 1,3}, TIMEX 1.4) 
1,TIME(1S), TIME( 1,6), TIME(1,7), TIME( 1,8), TIME( 1,9) 
1,TIMEQ, 10) 
RBAD (20,1 10)SPDY1, 1) SPD(L2),SPD(1,3).SPD¢ 1.4), SPD(1,5) 
1, SPD(1,6),SPD(1,7)SPDC_8),SPOX1,9)SPD(1, 10) 
READ (20, {20)R TIME{ 1.2), R TIMEX 1,2). R TIME( 1.3) TIME(1.4) 
1,RTIME(1,S)R TIME(1,6), RTIME( 1,7) 2 TIME(] 8). TIME( 1,9) 
1 RTIME(1,10} 
READ (20,130)GTIME(1, 1), GTIME(1,2),GTIME( 13), GTIME(1,4) 
1,GTIMEQ.5),GTIME(1,5),G TIME(], 7), TIMEX 1, 8), GTIME(1,9) 
LGTIME({1,10) 


OPEN (21, FILE = "P2LNK od”) 
READ (21,1 10/TIMEQ, 1} TIMEQ 2), TIMEG:3) TIME(,4) 

1, MEQ, 5), TIMB(2,6), TIMEQ, 7), TIME(2,8), TIME(2,9) 
 TIMEQ2, £0} 

READ @1, LIOSPD(2, 1) SPD(2,2),.SPD(2,3) SPDQ2,4),SPDQ.5) 
1,SPD(2,6),SPD(2, 7), SPD2,8),SPLX2, 9) SPD(2, 10) 

READ (2!,120)R TIMEX2, 1) R TIME(Z 2) 8 TIME(2,3) RTMELZ, 4) 
L.RTIMEQ, 5S) RTIMEQ, 6) RTIMEQ, 7) R TIME(2, 8) RTIME(2, 5) 

1 ROIME(, 10) 

READ (21,130)GTIMEQ, |), GTIME(2,2), GTIME(2,3), GTIMEX2,4) 
1,GTIME(2,5),GTIME(, 6),G TIME(, 7),G TIME(2, 8), GTIME(2,9) 
L.GTIMEX1, 10) 


OPEN (22, FILE = "PSLNK txt") 
READ Q22,! 10}TIME@, |), TIME(3, 2), TIMEG3,3) TIMEG,4) 
1, TIMEG, 5), TIMEG,6), TIMER, 7) TIMBG,8) TIME(S 9) 
1, TIMEQ, 10) 
READ 22,1 I)SPDG, 1), SPDG,2),SPDG,3),SPDG,4),SPDG,S) 
L.SPD(3.6)SFDG,7)SPDG.8),SPD(3,9,SPDG, 10) 
RBAD (22,1 2)R TIME(S, 1) RTIMEG, 22 TIMEG,3)R TIMER, 4) 
1. RTMEG, 5)RTIMEG 6) RTIME(3, 7) R TIMES §),RTIMEG3,9) 
LRTMEG, 10) 
READ 22, 130GTIMEG, 1), GTIME(, 2} GTIMEQ,3)GTIME3,4) 
1,GTIMEG,5), GIMEG,6),GTIME, 7.GTIMEG, 8), GTIME3.9) 
L.GTIMEQ,10) 

c 
OPEN (23, FILE = "P4LNK tx1”) 
READ (23, 10) TIME(4 1) TIME(4,2), TIME(43) TIME(4,4) 
1, TME(4,5), IMBG,6) TIMB(4, 7), TIME(43), TIME(4,9) 
1, TIME(4, 10) 
READ (23,110)SPD(4,1), SFD@,2).SPD(4,3)SPD(4,4),SPD(4.5) 
‘LSPDE4. 6), SPD(4,7), SPD, 8), SPD(4.9 SPPX, 10) 
READ (23,120) TIME(4 1) RTIME(4.2)R TIMEY43),8 TIME(4,4) 
1. RTIME(4,5)R TIME 6) R TIME(S, 7). TIME(4 8). R TIM 49) 
LRTME(,10) 
READ (23,1 30)GTIME(4, 1) GTIME(4,2),G TIME(4,3),GTIMEX44) 
L,GTEMEY4,5),GTIME(4 6} TEME(4, 7) GTIME(4 B), GTIME(4 9) 
LGTIME(4,10) 


OPEN (24, FILE = "PSLNK ut”) 
RBAD (24,1 1O)TIME(S, 1), TIME(S, 2), TIME(S 3), TIMES, 4) 
1, TIMES, 5), TIME(S, 5} TIMBG, 7}, TIMEXS 8), TIME(5,9) 
1TIMEG,10) 
READ (24, 1OSPD(S. 1). SPDS.2).SPDGS 3) SPD(5.4)SPOXS,5) 
1 SPDG,O,SPD(5,7), SPOS.) SPINS, I), SPDC5, 10) 
READ (24,1200 TIME(S, 1), A TIME(S 2) R TIME(S3). R TIME(S, 4) 
LRTIMEG, 5). TIMEG, 6). RIMES, 7). TIME(S, 8) RTIMELS,9) 
LRIIMEG, 10) 
READ (24,130)GTIMELS, 1), GTIME(S,2).GTIME(5,3),G TIMEX 5,4) 
1,GTIMESS 5), GTIMES5, 6), GTIMELS, 7), GTINEXS, 8) GIIME(S.9) 
LGTIMBGS, 10) 

110 FORMAT (10°73) 

120 FORMAT (100(419)) 

130 PORMAT (10(F14.3)) 

C SORT GPS TIME FROM LASTEST TO EARIEST 

OPEN (0, FILE="AVGLTT.CSV") 
DEX!) = (5843, <28/SZR0) 
DIZ) = (16725 5VS200) 
DD) = (1139. 12575280) 
DD) = (1115839659280) 
DIS) = (838.457/5280) 


DDO = 2176.219/5280) 
DD?) = (1134,646/5280) 
DDB) = (144216165220) 
DD) = (1257.546/5280) 
DD 10) = (6407 B6/S280) 
DOJ= 110 
TIMEQ)= 0.0 
NQ)=00 
LASTRUNG@)~ MAX(GTIME( 1,3) GTIME(2,)),.G TIME J). GTIME(4,J) 
1,GTIMEGS)) 
DOI=1,5 
¢ WITHIN 15 MINUTES(%00 SECONDS) AFTER LAST RUN 
FF (LASTRUNG)}GTIMEGJ) LA 900 sxxd. érne(ij) ne. 0) THEN 
THIMEQ) = TTIMED)} UME) 
‘TspdQ) = ep) + sok) 
NQ)=NOyH.0 
dae 
tienieg}) = tamol)) 
tepd) = tad) 
nf)=0g) 
END IF 
IF (LASTRUNQ)-GTIME(L) BQ. 0) THEN 
IROMED)“ ROME) 
END JF 
END DO 
€ CALCULATE AVG TIME MEAN SPEED OF VEHICLE WITHIN LSMENS. 
RANGE 
ATTIMEQ)= THIMEJYN) 
ASPEED(Q) — (DDG@))*60/A TTIMEQ) 
if Gettimea()) eq. 0.0) then 
speed) = 0.0 
adif 
if(e) 24.0.0) then 
mapeaxi{J) - 0.0 
wuime(j)= 0.0 
ond if 
if (o@)-1.0 EQ. 0) then 
ASPEEDG) tapi) 
endif 
END DO 
c 
WRITE(30, (A, 10(7,",A)) TIME(US40-US6)*, “US 40(0F FH SAQ(ON)? 
1,"US40(0N)-GRG TWN(OFFY" 
12GRG TWNOFF)GRG TWN(ON)" 
1"GRGB TWN(ON)SLVR PLM(OFFY 
1,°SLVER PLM(OFF}-SLVR PLMYON)” 
1,°SLVR PLM(ON)-BKRVLLOOFFY" 
1,"BKRVLL(OFP)-BERVLL{ON)* 
17SKRVLL{ON}-HRMN GLCHIOFFY” 
1,"ARMN GLCH(GFF)HRMN GLCH(ON)" 
1,"HRMIN GLCERON)-URS(OFF* 
WRITEQ30,200) A TTIMEC),A TTIME(2),A TTIMEG)A TIIME(4} A TIME(S) 
1 ATIIME(6)A TTIME(7).A TTIME(8).A TDIME(9),A TTIME( 10) 
WRITE(30,200) ASPHEDY( 1) ASPEED(2),ASPEED(3) ASPEEDX4) 
1, ASPEED(5),A SPEED(6),ASPEED(7),ASPEED(8) ASPEED(9),ASPEED( 10) 
WRITE30,200)N( NZ NGNONONONTNENAN(1 0) 
WRITE(30,201) IRTIME(1) 
WRITE(30,201)1R TIME(2) 
WRITE(0,201) IR TIMEG) 
WRITEGO,201)IR TIME(4) 
WRITE(30,20 IY ERTIMB(S) 
WRITE(30,201)1R TRME(S) 


WRITEG30,201)1R TIME(7) 
WRITE(0,201JTR TIME(S) 
WRITE(30,201)IRTIME() 
WRITEO201)IRTIME(!0) 

M0 FORMA T(I0C,"F7.3)) 

201 FORMAT (A,"AMPM”) 
END SUBROUTINE 
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APPENDIX D: 
SAMPLE GPS DATA 


Bob ileKey", obit eID", "RapecdT ime”, “PosT int ime", 9TTasec", Lavioude", "Longitude", "Ait itoda”, Speed”. 
course", ‘PopageFlag®, “¢psPixSrc". “Status”, *Seatrusz", “Seavusl ive", “MaminOil*, *Ond£2Z1". "Mn0f£2® , “Onb£t3 
") "OnOLEb", "OnQ£ES", "OnO£fh*, "Ond£E7", "On0SES", "Baneading™, "AeadingL"”, "Reading? *, "heading?" "Raadingt 
*, "LenExvended”, "Extended", Event ID", "Flagz” : 

¥78788,5 BAB LSLELS , SAB LE1320 , 0, 39754978 , “LV S1ALI7TS , LHOLLG ROTH, SSRI Lee Le Os eeccerscOeerecOr 20,0 
ITO TSED -£  C45LELG LS , 645151926 0, 39804039 ,~ LE LOSIEL, 2OGIL2 - 2290 GLI ibs ce“ AeDs reset rsDeres2Os 10/0 

278 VHO S$, CASAS LELS FES LS 1525.0 SIPBSEGS , - LOS LALIT, LEFB TO, LOSE BSI Lele ce LCs eect cee Ore rcv Oe 9.0 
JIBVEL, 2 SABLILELS - BSE LELIES 0299504087 ,-LOBLISAOY- 909200 -BBEL, OFZ Lil sre LOs cs sarersbaesreGe 28,0 
TISISL S$, FABLELSLD , C45 LELISO , O , FP 7OFI4A, - LOT LABOOS , LETZLS LOIS, ISPS Les ~ Le Oe esacvea Griese 202600 
JIBIES 2 - CQSLSLS LB , 645151390 ,0, 39804600 , - 106191584, 207004, 3120 693 Lele r ee LeQa aes anee Verse Be 08 
ITBIS4, 5, BAGLSLSLG , FAS LI LEW ,0, 97789894, - 109149097, LSSFET BEST OSVO Lek ee LB rec ce sr srerer Or e0c0 
PIBISS ,&, FHLSLSTS , 645151940 0, 39504861, ~-L0FLO7955, 200080 BUSA, FOZ, Lif, 2 be Oy es cnree Gere eeDe 20,9 
FIOWE 5 , BABLELSSO , CS LSL348 , 0, ISTIOMTL, - 105149054, SELL, 2620, 9598, Lk eke Orerer reer Orerced, 0,0 
2TET07 2 64151596 , 646151945 ,0,39504309 77208168152, 307917 3142 .1033,1-2,. 71,6, ocenecaBrerceO, 0,8 
TPB HS ,L  C4HLELSIO , SABLILESO GO -FHECAOCT , - LOS LOSIES ,, 206279 -FOWO, LOAF Lele se Le Or cre recaeOrase r0e 4060 

3M T69 5, BAELS 1590 646151986 ,0, 39772047 ,-- LO5LRIOGY, LBSTIO, 2OTO-BSH5 Lele ge “Lads rer ss reeVeers G7s0,0 
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Figure E-1. Average Link Travel Time. 
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Figure E-2. Average Link Speed. 
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Figure E-4. Link Speed for Vehicle No. 4. 
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Figure E-5. Instantaneous Vehicle Speed for Vehicle No. 4. 
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